Hybrid protoneutron stars with the MIT bag model 



O. E. Nicotra, M. Baldo, G. F. Burgio, and H.-J. Schulze 

Dipartimento di Fisica e Astronomia, Universita di Catania and 
INFN, Sezione di Catania, Via Santa Sofia 64, 95123 Catania, Italy 
(Dated: February 5, 2008) 

We study the hadron-quark phase transition in the interior of protoneutron stars. For the hadronic sector, 
we use a microscopic equation of state involving nucleons and hyperons derived within the finite-temperature 
Brueckner-Bethe-Goldstone many-body theory, with realistic two-body and three-body forces. For the descrip- 
tion of quark matter, we employ the MIT bag model both with a constant and a density-dependent bag parameter. 
We calculate the structure of protostars with the equation of state comprising both phases and find maximum 
masses below 1.6 solar masses. Metastable heavy hybrid protostars are not found. 
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I. INTRODUCTION 



^ , It IS generally believed that a neutron star (NS) is formed as a result of the gravitational collapse of a massive star (M ;> SM©) 
in a type-II supernova Just after the core bounce, a protoneutron star (PNS) is formed, a very hot and lepton-rich object, 

where neutrinos are temporarily trapped. The following evolution of the PNS is dominated by neutrino diffusion, which results 
first in deleptonization and subsequently in cooling. The star stabilizes at practically zero temperature, and no trapped neutrinos 
are left. 

, The dynamical transformation of a PNS into a NS could be strongly influenced by a phase transition in the central region of 
the star Calculations of PNS structure, based on a microscopic nucleonic equation of state (EOS), indicate that for the heaviest 

fSJ PNS, close to the maximum mass (about two solar masses), the central particle density reaches values larger than 1 /fur". In this 
density range the nucleon cores (dimension « 0.5 fm) start to touch each other, and it is likely that quark degrees of freedom 

OO will play a role. 

In a previous article f3\ we have studied static properties of PNS assuming that nucleons, hyperons, and leptons are present in 
stellar matter. Our calculations are based on the EOS derived within the Brueckner-Bethe-Goldstone (BBG) theory of nuclear 
matter, extended to finite temperature. We have found that for purely nucleonic stars, both thermal effects and neutrino trapping 
G slightly soften the EOS, thus reducing the value of the maximum mass. If hyperons are included, neutrino trapping shifts their 
onset to larger baryon density, and instead stiffens the EOS. This could lead to metastability during the subsequent evolution of 
Q the PNS to the late neutrino-free stage. 

^ In this work we extend the previous calculations, and take into account a possible hadron-quark phase transition. In fact, as 

c/2 in the case of cold NS, the addition of hyperons demands for the inclusion of quark degrees of freedom in order to obtain a 
maximum mass larger than the observational lower limit. For this purpose we use the BBG EOS for describing the hadronic 
J> phase and the MIT bag model at finite temperature for the quark matter (QM) phase. We employ both a constant and a density- 
dependent bag parameter B. We find that the presence of QM increases the value of the maximum mass of a PNS, and stabilizes 
it at about 1.5-1.6 Mq, no matter the value of the temperature. In contrast to purely hyperonic stars, neutrino trapping in hybrid 
stars does not increase the maximum mass and does not allow metastable states. 

The paper is organized as follows. In section|ll]we review the determination of the baryonic EOS comprising nucleons and 
hyperons in the finite-temperature Brueckner-Hartree-Fock approach. Section|ni]concerns the QM EOS according to the MIT 
bag model. In section CVl we present the results regarding PNS structure combining the baryonic and QM EOS for beta-stable 
nuclear matter. Section|V]contains our conclusions. 



II. BRUECKNER THEORY 

A. EOS of nuclear matter at finite temperature 

In the recent years, the BBG perturbative theory for nuclear matter has made much progress, since its convergence has 
been firmly established and it has been extended in a fully microscopic and self-consistent way to the hyperonic sector 
iill. Only few microscopic calculations of the nuclear EOS at finite temperature are so far available. The first semi- 
microscopic investigation of the finite-temperature EOS was performed in Ref. 1 8]. The results predict a Van der Waals behavior 
for symmetric nuclear matter, which leads to a liquid-gas phase transition with a critical temperature 7^ « 18-20 MeV. Later, 
Brueckner calculations L9> .101 and chiral perturbation theory at finite temperature 1 11 ] confirmed these findings with very similar 
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values of 7^. The Van der Waals behavior was also found in the finite-temperature relativistic Dirac-Brueckner calculations of 
ll2t.l3J . although at a lower temperature. 

The formalism which is closest to the BBG expansion, and actually reduces to it in the zero-temperature limit, is the one 
formulated in In this approach the essential ingredient is the two-body scattering matrix K, which, along with the single- 
particle potential U, satisfies the self-consistent equations 

{k,k2\K{W)\k3k4) = {kik2\V\k3k4) 

+ Re ^ {k, k2 \V m) ^w_"l'^^^iV"rH^ {kmW)^) 



,77, W -E,,, - Eu + ie 



(1) 



and 



U{ki)=Y,n{k2){kik2\K{W)\kik2)A, (2) 

where ki generally denote momentum, spin, and isospin. Here V is the two-body interaction, and we choose the Argonne Vig 
nucleon-nucleon potential 1 1 5 ] . W =Ei^^ + Ej^^ represents the starting energy, — k^j 2m + U{k) the single -particle energy, and 
n{k) is a Fermi distribution. At T = Eq. Q coincides with the Brueckner equation for the K matrix at zero temperature. 

For given nucleon densities and temperature, Eqs. Q and (|3 have to be solved self-consistently along with the equations for 
the densities, p; — LA:"i(^)' the free energy density, which has the following simplified expression 



/= L 



k 



(3) 



where 



Si = -Y,(niik)\nn,ik) + [1 -«,(A:)]ln[l (4) 

is the entropy density for component / treated as a free gas with spectrum Ei{k). 

In deriving Eq. (|3}, we have introduced the so-called Frozen Correlations Approximation, i.e., the correlations at T 7^ are 
assumed to be essentially the same as at T = 0. This means that the single-particle potential Ui{k) for the component / can be 
approximated by the one calculated at T = 0. This allows to save computational time and simplify the numerical procedure. It 
turns out that the assumed independence is valid to a good accuracy, at least for not too high temperature (see Ref. Id Oil . Fig. 12). 
For a more extensive discussion of this topic, the reader is referred to Ref. 1 16], and references therein. 

In our many -body approach, we have also introduced three-body forces (TBF) among nucleons, in order to reproduce correctly 
the nuclear matter saturation point po ~ 0.17 fm^^, ii/A « — 16 MeV and obtain values of incompressibility and symmetry 
energy at saturation point compatible with those extracted from phenomenology ||T7ll . Since a complete microscopic theory of 
TBF is not available yet, we have adopted the phenomenological Urbana model 118], which in the BBG approach is reduced to 
a density-dependent two-body force by averaging over the position of the thkd particle, assuming that the probability of having 
two particles at a given distance is modified according to the two-body correlation function fulfl^ . 

The fast rise of the nucleon chemical potentials with density in NS cores f2^,'2l'l may trigger the appearance of strange 
baryonic species, i.e., hyperons. For this purpose we have extended the BBG approach in order to include the and A hyperons 
in a fully self-consistent way |5l01- In our work we have used the Nijmegen soft-core nucleon-hyperon (NH) potentials NSC89 
ll22il . and neglected the hyperon-hyperon (HH) interactions, since so far no reliable HH potentials are available. Recently also 
calculations with the NSC97 potentials were concluded |7|, which do include extensions to the HH sector, but yield very similar 
results. 

The presence of hyperons strongly softens the EOS, mainly due to the larger number of baryonic degrees of freedom. This EOS 
produces a maximum NS mass that lies slightly below the canonical value of 1.44 Mq |23], which could indicate the presence 
of non-baryonic (quark) matter in the interior of heavy NS [24, 25, 26]. However, the quantitative effects of more reliable NH 
and HH potentials as well as of hyperonic TBF need still to be explored in the future. For these reasons, we will in this article 
perform finite-temperature calculations in a much simpler way by using non-interacting hyperons. This approximation is well 
justified for the present work, because the hadron-quark phase transition occurs at very low density, where hyperons do not yet 
play a role. 

In any case, our microscopic baryonic EOS turns out to be very soft and features very different characteristics in comparison 
to the often used relativistic mean field (RMF) models. In particular, the combination of a QM phase with a RMF model leads 
usually to a reduction of the maximum NS mass Emu, whereas in our approach the maximum mass increases due to the 
presence of QM, which is in fact required in order to cover the current observational values of NS masses. Also for PNS very 
different properties will arise, as we will see. 
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Baryon density pg (fm ') 

FIG. 1: The EOS in beta-stable baryon matter is displayed for temperatures T = and 30 MeV. The upper (lower) panel displays results for 
neutrino-free (neutrino-trapped) matter, for purely nucleonic (upper curves) and hyperonic (lower curves) stellar matter. 



B. Composition and EOS of hot stellar matter 

For the determination of the composition of beta-stable baryonic matter, one needs to know the functional dependence of the 
free energy, Eq. (|3}, on the individual partial densities and on temperature. In Ref. |3] we have provided analytical parametriza- 
tions of our numerical results for this purpose. From the free energy density one can then calculate the nucleon chemical 
potentials 

df 

= ^ , (5) 

whereas the chemical potentials of the noninteracting leptons and hyperons are obtained by solving numerically the free Fermi 
gas model at finite temperature. 

For stars in which the strongly interacting particles are only baryons, the composition at given baryon density p = ^, p, is 
determined by the requirements of beta equilibrium and charge neutrality, i.e., 

I^i = biH„-qi{iii-iJ.y,), (6) 

= J^qiXi + Y,qiXi, (7) 
i I 

where x, = p,7p is the baryon concentration, bj the baryon number, and qi the electric charge of the species /. Equivalent 
quantities are defined for the leptons I — e,jj.. The initial PNS contains trapped neutrinos, so the electron and muon lepton 
numbers 

Yi =xi-xi+xvi ~Xvi (8) 

are conserved on dynamical time scales and we fix them to = 0.4 and — 0, as indicated by gravitational collapse calculations 
of the iron core of massive stars. 

IP-* 

Solving these equations, we have found in |3] that the electron fraction is larger in neutrino-trapped than in neutrino-free 
matter, and, as a consequence, the proton population is larger. Moreover, neutrino trapping shifts the threshold density of the 
to high density, whereas A's appear at slightly smaller density. This is due to the fact that the onset depends on the 
neutron and lepton chemical potentials, i.e., fii^- = jj.,, + jJ-e — piv^, which stays at larger values in neutrino-trapped matter than 
in the neutrino-free case because of the larger fraction of electrons, thus delaying the appearance of the E^ to higher baryon 
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density and limiting its population to a few percent. On the other hand, the A onset depends on the neutron chemical potential 
only, /Xa = ix„, which stays at slightly lower values in the neutrino-trapped case. When the temperature increases, the hyperon 
thresholds disappear and more and more hyperons are present also at low densities, but they represent only a small fraction of 
the total baryon density in this region of the PNS. Altogether, the hyperon fractions are much smaller than in the neutrino-free 
matter and therefore the corresponding EOS will be stiffer than in the neutrino-free case. 

Once the composition of the beta-stable stellar matter is known, one can proceed to calculate the free energy density / and 
then the pressure p through the usual thermodynamical relation 

P = P'^. (9) 

The resulting EOS is displayed in Fig.Q] where the pressure for beta-stable stellar matter, without (upper panel) and with (lower 
panel) neutrinos, is plotted as a function of the baryon density at temperatures T — Q and 30 MeV. We notice that thermal effects 
produce a slightly stiffer EOS with respect to the cold case, and that at very high densities they almost play no role. The inclusion 
of hyperons, however, causes a dramatic effect, because the EOS gets much softer, no matter the value of the temperature. In 
the neutrino-trapped case, the EOS is slightly softer than in neutrino-free matter if only nucleons and leptons are present in the 
stellar matter. Also here the presence of hyperons introduces a strong softening of the EOS, but less than in the neutrino-free 
case, because now the hyperons appear later in the matter and their concentration is lower This fact could lead to metastable 
stars which suffer a delayed collapse while cooling down, as discussed in fzfl . 



III. QUARK MATTER EQUATION OF STATE 
A. The MIT Bag Model 

We review briefly the description of the bulk properties of uniform QM at finite temperature, deconfined from the beta-stable 
hadronic matter discussed in the previous section, by using the MIT bag model [2m- In its simplest form, the quarks are 
considered to be free inside a bag and the thermodynamic properties are derived from the Fermi gas model, where the quark 
q = u,d,s baryon density, the energy density, and the free energy density are given by 

^ 1/(03 [/.^« ^10) 

= sLJ [/.^W W] E,{k)+B, (11) 
/e = £Q-TY.Sg, (12) 



where g = 6 is the quark degeneracy, Etj{k) = y + k^, B is the bag constant, Sq the entropy density of a noninteracting quark 
gas, and the Fermi distribution functions for the quarks and anti-quarks are 

l+exp[{E,{k)TH,)/T] ^^^^ 

with jig being the quark chemical potentials. The corresponding expressions at T = can be obtained by eliminating the 
antiparticles and substituting the particle distribution functions by the usual step functions. We have used massless u and d 
quarks, and iris = 150 MeV. 

It has been found that within the MIT bag model (without color superconductivity) with a density-independent bag 

constant B, the maximum mass of a NS cannot exceed a value of about 1.6 solar masses. Indeed, the maximum mass increases 
as the value of B decreases, but too small values of B are incompatible with a hadron-quark transition density p > 2-3 po in 
nearly symmetric nuclear matter, as demanded by heavy-ion collision phenomenology. Values of B > 150 MeV/fm^^ can also 
be excluded within our model, since we do not obtain any more a phase transition in beta-stable matter in combination with our 
baryonicEOS |24]. 

In order to overcome these restrictions of the model, one can introduce a density-dependent bag parameter B{p), and this 
approach was followed in Ref. i24ll . This allows one to lower the value of B at large density, providing a stiffer QM EOS and 
increasing the value of the maximum mass, while at the same time still fulfilling the condition of no phase transition below 
p w 3po in symmetric matter. In the following we present results based on the MIT model using both a constant value of the bag 
parameter, B — 90 MeV/fm^, and a gaussian parametrization for the density dependence. 



B(p) =Boc + (Bo-fioc)exp 



Vpo 



(14) 
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FIG. 2: Particle fractions in beta-stable QM at T = (left panels) and T = 50 MeV (right panels) for neutrino-free (upper panels) and 
neutrino-trapped (lower panels) matter 



with Boo = 50MeV/fm-3, Bq = 400MeV/fm3, and j3 = 0.17, see Ref. UM- 

The introduction of a density-dependent bag has to be taken into account properly for the computation of various thermody- 
namical quantities; in particular the quark chemical potentials and the pressure are modified as 

dp 

Nevertheless, due to a cancelation of the second term in (I15> . occurring in relations illi for the beta-equilibrium, the composition 
at a given total baryon density remains unaffected by this term (and is in fact independent of B). At this stage of investigation, 
we disregard possible dependencies of the bag parameter on temperature and individual quark densities. For a more extensive 
discussion of this topic, the reader is referred to Refs. J^HEES- 



B. Quark matter in beta equilibrium 

In a PNS with QM and trapped neutrinos we must add the contribution of the leptons as free Fermi gases to the (free) energy 
density, Eqs. ( II lll2t . and impose beta equilibrium, charge neutrality, and baryon and lepton number conservation I2lll . 
More precisely, the individual quark chemical potentials are fixed by Eq. (|6} with bg =1/3, which implies 

11,1 = H, = H„ + Hi-livi. (17) 

The charge neutrality condition and the total baryon number conservation read 

= 2p^-p,i-p,-pe, (18) 

P = Pu + Pcl + Ps, (19) 

and Eq. specifies lepton number conservation. These equations determine the composition Pq{p) and the pressure of the QM 
phase, 

PQ{p)=P^-fQ. (20) 
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FIG. 3: Baryon density (lower panels) and pressure (upper panels) as a function of baryon chemical potential of beta-stable baryonic matter 
(solid curves) and quark matter (dashed curves) with (right panels) and without (left panels) neutrino trapping at different temperatures T = 
0, 10,30,50 MeV. The vertical dotted lines indicate the positions of the phase transitions. A bag constant B = 90 MeV/fm-' is used for QM. 




FIG. 4: Same as Fig.|3| but with a density-dependent bag parameter. 



In Fig.|2lwe plot the particle fractions jc, = p,/ p as a function of baryon density for neutrino-free (upper panels) and neutrino- 
trapped (lower panels) QM. Since in the range of temperature considered here thermal effects are rather weak, we report only 
results for T = and T =50 MeV. Conversely, the presence of neutrinos influences quite strongly the composition: In this case 
the relative fraction of u quarks increases substantially from 33% to about 42%, compensating the charge of the electrons that 
are present at an average percentage of 25% throughout the considered range of baryon density, whereas d and s quark fractions 
are slightly lowered. One also can argue from the figure that most of the electrons present in the neutrino-free case come from 
thermally excited pairs. 

Before studying the phase transition inside PNS, we would like to reiterate the fact that our combination of baryonic and 
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FIG. 5: Free energy per baryon as a function of baryon density of beta-stable baryonic matter (solid lines) and quark matter (dashed lines) 
with a bag constant B = 90 MeV/fm-* (left panels) or a density-dependent bag parameter (right panels) with (lower panels) and without (upper 
panels) neutrino trapping at different temperatures T = 0, 10, 30, 50 MeV (upper to lower curves). 



QM EOS does not allow any phase transition in symmetric nuclear matter below about 3-4 po, as ruled out by observational 
evidence. This has been amply discussed in Refs. 13 1312^ and was actually an important motivation for our choice of the 
QM EOS, as mentioned above. 



C. Phase transition in hot beta-stable matter 



We now consider the hadron-quark phase transition in beta-stable matter at finite temperature. In the present work we adopt 
the simple Maxwell construction for the phase transition from the plot of pressure versus chemical potential. The more general 
Glendenning (Gibbs) construction i2lll30ll is still affected by many theoretical uncertainties lIsTIl and in any case influences very 
little the final mass-radius relations of massive (proto)neutron stars ll24ll . 

We therefore display in Figs. |3] and 0] the pressure p (upper panels) and baryon density p (lower panels) as functions of the 
baryon chemical potential /x„ for both baryonic and QM phases at temperatures 7" = 0, 10, 30, 50 MeV. The crossing points of the 
baryon and quark pressure curves (marked with a dot) represent the transitions between baryon and QM phases. The projections 
of these points (dotted lines) on the baryon and quark density curves in the lower panels indicate the corresponding transition 
densities from low-density baryonic matter, pn, to high-density QM, pg. The results in Fig.|3]are obtained with B = 90 MeV/ fm^ 
and those in Fig.0]with B{p). 

The following general observations can be made: 

(i) The transition density pn is rather low, of the order of the nuclear matter saturation density; 

(ii) The phase transition density jump pg — pn is large, several times p//; 

(iii) Thermal effects shift pn to lower values of subnuclear densities and increase the density jump pg — pn', 

(iv) Neutrino trapping lowers even more the transition density. For the cold case the presence of neutrinos even inhibits com- 
pletely the phase transition. However, this feature has no physical relevance, because no trapping occurs in cold catalyzed 
neutron stars. The reason for such behavior is the stronger increase of the QM pressure due to neutrino trapping compared to the 
one of the baryonic phase, as evidenced in the figure. This leads naturally to an earlier onset of the QM phase; 

(v) The model with density-dependent bag parameter predicts (by construction, see the discussion above) larger transition den- 
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2 4 6 8 10 20 40 60 80 100 
P / Po R (km) 

FIG. 6: (Proto)Neutron star mass-central density (left panel) and mass-radius (right panel) relations for different core temperatures T = 
0, 10, 30, 40, 50 MeV and neutrino-free (solid curves) or neutrino-trapped (dashed curves) matter A bag constant fi = 90 MeV/ fm-' is used for 
QM. 
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FIG. 7: Same as Fig.|6| but with a density-dependent bag parameter. 

sities Ph and larger jumps pg — than those with bag constant B = 90 MeV/fm^. The plateaus in the Maxwell construction 
are thus wider for the former case. 

These qualitative conclusions are reaffirmed in Fig.|5j where we report the free energy perbaryon, F/A, as a function of baryon 
density, the four panels corresponding to the cases mentioned above with and without neutrino trapping and with different bag 
parameters. Here the crossing points between solid lines (baryon phase) and dashed lines (QM phase) indicate an average phase 
transition density, above which the QM phase is characterized by a value of the free energy lower than that of the baryon phase, 
thus indicating that the QM phase is energetically favoured. 

Based on these results for the beta-stable baryon and QM phases, we proceed now to the determination of the properties of 
static (proto)neutron stars. 

IV. PROTONEUTRON STAR STRUCTURE 

We assume that a (proto)neutron star is a spherically symmetric distribution of mass in hydrostatic equiUbrium. The equilib- 
rium configurations are obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) equations fl'l for the pressure p and the 
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FIG. 8: Pressure as a function of energy density for beta-stable matter witli (dashed curves) and without (solid curves) neutrino trapping at 
different temperatures T = 0, 30, and 50 MeV with a bag constant S = 90 MeV/fm-' (left panel) or a density-dependent bag parameter (right 
panel). 




FIG. 9: Relation between gravitational mass and baryonic mass for cold neutron stars (solid curves) and protoneutron stars at T = 30,40 MeV 
(dashed lines) with a bag constant fi = 90 MeV/fm-' (upper panel) or a density-dependent bag parameter (lower panel). 



enclosed gravitational mass m and baryonic mass mg, 

dp{r) Gm{r)e{r) 



dr 



Jl+p{r)/e{r)] [l+47rrV(r)/m(r)] ^^^^ 
1 ^2Gm{r)/r 

ATtr^e{r) , (22) 



dm{r) 
dr 

dniB [r) Anr^p (r) 



dr y/l-2Gm{r)/r' 



(23) 



where niN = 1 .67 x lO^^^g is the nucleon mass and G the gravitational constant. Starting with a central mass density e(r = 0) = 
Ec, we integrate out until the pressure on the surface equals the one corresponding to the density of iron. This gives the stellar 
radius R and the gravitational mass Mq = m{R) and baryonic mass Mb = niBiR). 

Moreover, in our model we assume that a PNS, in its early stage, is composed of a hot isothermal and neutrino-opaque core 
separated from an outer cold crust (described by the Baym-Pethick-Sutherland |32] and the Feynman-Metropolis-Teller i33ll 
EOS) by an isentropic, beta-equilibrated, and neutrino-free intermediate layer in the range of baryon density from 0.01 fm^^ 
down to 10^^ fm^^ |34], which is based on the EOS LS220 of Lattimer and Swesty [S^J with compressibility K — 220 MeV. 
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For the isothermal core we use as input the BHF EOS for the baryonic matter and the MIT bag models for the beta-stable QM 
phase, as discussed above. In order to ensure an exact matching of all the thermodynamic quantities between core and crust we 
perform a fine tuning of the entropy per baryon of the isentropic envelope, 5env., around the matching value of baryon density 
Penv. ~ 0.01 fm^^. More precisely, fixing San. to the values 6,8, 10 (in units of the Boltzmann constant) we get temperature 
profiles which drop quickly from Tcoie = 30,40,50 MeV to zero in the considered range of density, thus suggesting a natural 
correspondence between values of Tcore and Sem. 1 34]. 

In other words, we model the PNS by adopting a hybrid isothermal plus isentropic temperature profile, allowing us to take into 
account as much as possible information coming from dynamical simulations |27, 36, 37, 38], which show that the temperature 
drops rapidly to zero at the surface of the star due to the fast cooling of the outer part of the PNS, where the stellar matter is 
transparent to neutrinos. In addition, during the early stage, the outer part of a PNS is characterized by a high value of the 
entropy per baryon. Essentially, in the first ten seconds, the entropy profile decreases from the surface to the core starting from 
values of 6-10 |36]. 

We consider this approach more realistic than the one employed for baryonic stars in |3], where a temperature profile in the 
shape of a step function was assumed, joining directly the baryonic EOS and the one for the cold outer crust. This procedure 
leads, however, to unstable PNS at Tcore slightly above 30 MeV and also yields a strong dependence on the value of density 
where the cold crust is attached of the minimum PNS mass, which is a fundamental quantity characterizing the range of stability 
of such hot and compact objects. The improved procedure described above applied to baryonic PNS leads to more reliable values 
of the minimum gravitational mass, while hardly affecting the value of the maximum mass, which mainly depends on the EOS 
employed and the core composition. For the purely nucleonic case (no hyperons) we find now a maximum mass within a narrow 
range of 1.77-1.87 l34ll. similar to jl. 

In the following we present results schematizing the entire evolution of the star as divided in two main stages. The first consists 
of a PNS with a hot (Tcore ~ 30-50MeV) neutrino-trapped core and a high-entropy transition layer (Senv. ~ 6-10), joined to a 
cold outer crust. The second stage represents the short-term cooling, where the neutrino-free core possesses a low temperature 
of about 10 MeV and is direcly attached to a cold crust at p = 3 x 10^^ fm^^. 

A. Numerical Results 

The results are plotted in Figs. |6l and where we display the gravitational mass Mq (in units of the solar mass M0) as a 
function of the radius R (right panels) and the central baryon density Pc (left panels), for QM EOS with B = 90 MeV/fm^ and 
Z?(p), respectively. 

Due to the use of the Maxwell construction, the curves are not continuous f27]: For small enough central densities (large 
radii) the stars are purely baryonic. Then a sudden increase of the central density is required in order to initiate the QM phase 
in the center of the star, corresponding to the phase diagrams Figs.|3]and|3 By performing the Glendenning construction, the 
curves would become continuous. Heavy PNS in our approach are thus practically quark stars with only a thin outer layer of 
baryonic matter. Below the maximum mass configuration, however, the stars develop an extended outer envelope of hot matter, 
the details of which depend on the treatment of the low-density baryonic phase and the phase transition. We will therefore in the 
following focus on the properties of heavy stars close to the limiting mass, which are unaffected by these complications. 

For completeness we display the complete set of results at core temperatures T = 0, 10,30,40,50 MeV with and without 
neutrino trapping, although only the curves with high temperatures and neutrino trapping and low temperatures without trapping 
are the physically relevant ones. We observe in any case a surprising insensitivity of the results to the presence of neutrinos, in 
particular for the = 90 MeV/fm^ case, which can be traced back to the fact that the QM EOS p{e) is practically insensitive 
to the neutrino fraction. In fact, for the B = 90 MeV/fm^ case and assuming massless quarks and leptons, the universal relation 
p = {e — 4B)/3 would hold irrespective of the internal composition of the quark-lepton phase. This is illustrated in Fig. |8] 
showing the beta-stable EOS p{e) at three temperatures with and without neutrinos. Indeed the pressure in the QM phase is 
nearly the same in both cases,even to a lesser degree for the B{p) EOS. 

On the other hand, the temperature dependence of the curves is quite pronounced for intermediate and low-mass stars, showing 
a strong increase of the minimum mass with temperature, whereas the maximum mass remains practically constant under all 
possible circumstances. Above core temperatures of about 40-50 MeV all stellar configurations become unstable. 

Concerning the dependence on the QM EOS, we observe again only a slight variation of the maximum PNS masses between 
1.55 Mq for B = 90 MeV/fm^ and 1.48 Mq for B{p). Clearer differences exist for the radii, which for the same mass and 
temperature are larger for the B = 90 MeV/fm^ model, as has also been found in 1,24. 26] for cold NS. The maximum PNS 
masses turn out to be very close to those of cold NS, thus excluding the possibility of metastable configurations. 

In order to illustrate better this issue, we show in Fig.|9lthe relation between baryonic and gravitational mass for cold NS and 
PNS at r = 30 and 40 MeV. In this plot the evolution of an isolated star proceeds on vertical lines connecting the upper curve for 
PNS with the lower curve for NS. One notes immediately that metastable hybrid stars do not exist in our model: For the B{p) 
case, the heaviest PNS (Mb = I.28M0, = I.48M0) transits into a Mq = \26Mq NS, which is consequently the heaviest 
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NS that can be produced without accretion in our approach. Accretion might then further augment the NS mass to a maximum 
stable value of Mq ~ \.53Mq,. 



V. CONCLUSIONS 

In this article we extended our previous works on cold baryonic f6'l and hybrid fl^,^2?^2^ NS and baryonic PNS to the 
case of hybrid PNS. We combined the most recent microscopic baryonic EOS in the BHF approach involving nuclear three- 
body forces and hyperons with two versions of a generalized MIT bag model describing the QM phase: One using a fixed bag 
constant, the other one a density-dependent bag parameter B(p ) in order to explore the maximum PNS mass that can be reached 
in this approach. 

We modelled the profile of a PNS in an extremely simplified way with a constant lepton fraction Yg = 0.4 and constant temper- 
ature of the core and an isentropic cover, leaving the core temperature as a global parameter. A really satisfying treatment would 
require coupled dynamical simulations for the various microscopic and macroscopic evolution equations, which is currently 
beyond our reach. 

We found in Ref. 1 3] that purely baryonic (hyperonic) PNS can reach masses of about 1 .5 Mq, and nearly the same mass limit 
is now obtained also for hybrid PNS. The difference between both configurations lies in the transition to a cold NS, which for 
baryonic stars could lead to a delayed collapse to a black hole, as the cold baryonic NS cannot support masses up to 1 .5 Mq 0], 
whereas this phenomenon is exluded for hybrid NS, which can sustain masses up to those of the PNS. 

This result is in contrast to Ref. |39|, where such metastability was found also for hybrid stars; however, with a much stiffer 
baryonic EOS and correspondingly larger value of the bag constant B = 200 MeV/fm^, which is excluded in our model. As 
mentioned before, our baryonic EOS is especially soft due to the presence of hyperons, compensating the repulsive character 
of nucleonic TBF at high density. Its associated maximum NS mass remains below 1 .4 solar masses, and the presence of QM 
inside the star is required in order to reach larger maximum masses. We have thus shown that the presumed increase of the 
hadron-quark phase transition density due to the presence of neutrino trapping |39, 40, 41] is not a general feature, but might be 
reversed in combination with a very soft baryonic EOS. 

Altogether, we once again confirm our prediction of rather low limiting masses for (proto)neutron stars, irrespective of varia- 
tions of the baryonic or QM EOS. Therefore, the experimental observation of a very heavy (M > 1 .8M0) NS, as claimed recently 
by some groups |42] (M « 2.2Mq), if confirmed, could hint to serious problems for the current rather simple theoretical mod- 
elling of the high-density quark matter phase, the assumptions about the phase transition between baryon and QM phase, or even 
for the TOV equations at extreme baryon densities. 
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